The UV Properties of Star Forming Galaxies I: HST 
WFC3 Observations of Very-high Redshift Galaxies 



Stephen M. Wilkins 1 *, Andrew J. Bunker 1 , Elizabeth Stanway 2 ' 3 , Silvio Lorenzoni 1 
Joseph Caruana 1 

1 Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, OX1 3RH, U.K. 

2 Department of Physics, University of Warwick, Gibbet Hill Road, Coventry, CV4 7AL, U.K. 

3 H H Wills Physics Laboratory, Tyndall Avenue, Bristol, BS8 1TL, U.K. 



1 July 2011 



ABSTRACT 

The acquisition of deep Near-IR imaging with Wide Field Camera 3 on the Hubble 
Space Telescope has provided the opportunity to study the very-high redshift Universe. 
For galaxies up to z « 7.7 sufficient wavelength coverage exists to probe the rest-frame 
ultraviolet (UV) continuum without contamination from either Lyman-a emission or 
the Lyman-a break. In this work we use Near-IR imaging to measure the rest-frame UV 
continuum colours of galaxies at 4.7 < z < 7.7. We are have carefully defined a colour- 
colour selection to minimise any inherent bias in the measured UV continuum slope 
for the drop-out samples. For the highest-redshift sample (6.7 < z < 7.7), selected as 
£/850Zp-band dropouts, we find mean UV continuum colours approximately equal to 
zero (AB), consistent with a dust-free, solar metallicity, star forming population (or 
a moderately dusty population of low metallicity) . At lower- redshift we find that the 
mean UV continuum colours of galaxies (over the same luminosity range) are redder, 
and that galaxies with higher luminosities are also slightly redder on average. One 
interpretation of this is that lower-redshift and more luminous galaxies are dustier, 
however this interpretation is complicated by the effects of the star formation history 
and metallicity and potentially the initial mass function on the UV continuum colours. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high-redshift - ultraviolet: galaxies 



1 INTRODUCTION 

The addition of Wide Field Camera 3 (WFC3) to the Hub- 
ble Space Telescope (HST) provides the sensitivity to effi- 
ciently probe the Universe at observed near-Infrared (NIR) 
wavelengths. These abilities are particularly important for 
studies of the very-high redshift Universe, as the rest-frame 
UV and optical emission is redshifted into the NIR. 

Broadband NIR imaging such as that obtained by 
WFC3 can be exploited by taking advantage of the red- 
shifted Lyman-a break (e.g. at z = 7.2 Lyman-a b s ~ 1/xm) 
to allow the identification of numerous high-redshift (6 < 
z < 9) star forming galaxy candidates (e.g. Oesch et al. 
2010, Bouwens et al. 2010b, Bunker et al. 2010, McLure et al. 
2010, Wilkins et al. 2010, Finkelstein et al. 2010, Wilkins et 
al. 2011, Bouwens et al. 2010c, Lorenzoni et al. 2011, Mclure 
et al. 2011). From these observations it is possible to con- 
strain the rest-frame ultraviolet (UV) luminosity function, 
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and thus the UV luminosity density. Making assumptions 
regarding the composition and properties of the stellar pop- 
ulations (i.e. metallicity, star formation history, initial mass 
function) and dust attenuation we can constrain both the 
star formation rate and ionising photon luminosity density 
(e.g. Wilkins et al. 2011). 

The deep WFC3 NIR imaging also allows the mea- 
surement of the rest-frame UV continuum colours of high- 
redshift 4.7 < z < 7.7 objects without the possibility of 
contamination from either the Lyman-a break or Lyman- 
a emission. The UV continuum colour is theoretically pre- 
dicted and empirically confirmed to be sensitive to the 
amount of dust reddening (e.g. Meurer et al. 1999, Pan- 
nella et al. 2009), potentially making it a useful diagnostic 
of the dust attenuation and thus providing a method for con- 
straining intrinsic star formation rate density (as opposed 
to the observed star formation density). Indeed several re- 
cent studies have investigated the UV continuum slopes of 
high-redshift galaxies (e.g. Stanway, McMahon & Bunker 
2005, Bouwens et al. 2009, Bunker et al. 2010, Bouwens et 
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al. 2010a, Dunlop et al. 2011) typically finding blue UV con- 
tinuum slopes consistent with negligible dust attenuation. 

Complicating this however is the fact that the UV con- 
tinuum is also affected by other properties intrinsic to the 
stellar population, including the previous star formation his- 
tory, metallicity, initial mass function and the presence of 
Lyman-a emission (unless it is outside the range of the UV 
continuum probed). While the sensitivity of the UV contin- 
uum colour to different influences complicates its interpre- 
tation as an estimate of dust attentuation, it may provide 
a technique to probe a range of other physical properties, 
especially at high redshift (e.g. Bouwens et al. 2010a). 

High-redshift galaxies are typically selected using a 
form of the Lyman-break technique (e.g. Steidel, Pettini and 
Hamilton 1995, Steidel et al. 1999) in which galaxies are 
principally selected on the basis of an extremely red colour, 
the cause of which is interpreted as the result of the Lyman-a 
break. This is typically combined with an additional colour 
selection and morphological criteria to guard against intrin- 
sically red contaminants, such as low-mass dwarf stars and 
lower-redshift evolved stellar populations (e.g. Stanway et 
al. 2008). However this makes studying the distribution of 
UV continuum colours difficult as most selection criteria de- 
signed to maximise the number of high-redshift candidates 
introduces both explicit and implicit biases. In this work we 
define a selection criterion which minimises these biases and 
use it to identify a catalogue of high-redshift star forming 
galaxies. Using this catalogue we investigate the rest-frame 
UV continuum colour distribution as a function of both lu- 
minosity and redshift. This paper is organised as follows: 
in Section 2 we describe the observations, data reduction 
and photometry. This is followed in Section 3 by a descrip- 
tion of our selection criteria and how it was derived, to- 
gether with simulations used to probe its effectiveness. In 
Section 4 we outline the various properties that affect the 
intrinsic UV continuum slope, and investigate how changing 
these properties (including the initial mass function and the 
star formation history) quantitatively affects the UV con- 
tinuum colour. In Section 4 we present our measurements 
of the UV continuum colour, together with the colour ex- 
cess E(B — V) (inferred assuming the Calzetti et al. 2000 
extinction curve) and the slope of the UV continuum /3 (i.e. 
fx oc A' 3 ) of our sample split by redshift and also luminosity. 
Finally in Section 5 we present our conclusions. Through- 
out this work we employ the AB magnitude system (Oke & 
Gunn 1983) and a ACDM cosmology (f2 A = 0.7, fi M = 0.3, 
and Ho = 70kms _1 Mpc -1 ) 



2 OBSERVATIONS, CATALOGUE 

CONSTRUCTION AND PHOTOMETRY 

In this work WFC3 observations of three fields (HUDF, 
P34 and P12, see Wilkins et al. 2011) in the vicinity of 
the CDFS/GOODS-South are combined with existing Ad- 
vanced Camera for Surveys (ACS) (6435™ (unavailable for 
the P34 and P12 fields), v 60 6w, hizw and z 850 i p ) imaging. 
The WFC3 data come from the HST Treasury programme 
GO-11563 (P.I. G. Illingworth) and consist of observations in 
three filters: Yiosu,, Ji25u> & Hieow A more detailed account 
of the observations, data reduction, and photometry is pre- 
sented in Wilkins et al. (2011) and Bunker et al. (2010). In 



each of the three fields the 2a point source detection limit 
is > 29.50 (AB) for the ACS images and > 29.15 for the 
WFC3 images with the exception of the Hisam imaging of 
the P12 field (~ 28.2). 



3 CANDIDATE SELECTION 

The combination of the HST ACS and WFC3 imaging avail- 
able allows us to identify star forming galaxies, and measure 
their rest-frame UV continuum colours from z = 2 — ► 8 (for 
the HUDF field). However, the addition of the new deep 
WFC3 NIR Yio5w, J125™, and Hieow imaging is most rele- 
vant to the highest redshifts (z > 4.5) and it is these objects 
on which we focus our attention. 



3.1 Lyman-Break Technique 

At sufficiently high redshift an increasing number of absorp- 
tion lines due to intervening neutral hydrogen causes a de- 
crease in the flux shortward of the rest-frame location of 
Lyman-a. At the lower redshift end this is manifested by 
the Lyman-a forest (Songaila & Cowie 2002); at very- high 
redshift the large number of absorption lines effectively effec- 
tively extinguishes the flux shortward of 1216A, producing 
the Gunn- Peterson trough feature (Gunn & Peterson 1965). 

In either case, the flux measured through broadband 
photometry shortward of rest-frame Lyman-a will be re- 
duced relative to that above (longward of) the break. This 
can be exploited by using at least two filters are positioned 
above and below the observed frame location of the Lyman- 
a break. 

Figure Q] shows the (z/ssoip — Y/ios™)ab colour of a 
galaxy with a simple, yet physically motivated spectrum 
(fx = A" at A rest > 1216A with 6 = -2.0 Q fx = 
0.01 x A" 2 at 912A < A re st < 1216A (i.e. intervening 
neutral Hydrogen absorbs ~ 99% of the emitted photons) 
and fx = 0.0 at A rest < 912A) as a function of red- 
shift. At z > 6. the Zf 85ozp filter encompasses the Lyman- 
a break resulting in a decrease in the recorded flux and 
an increasingly red (z f8 soi P - Y fl05w )AB colour. Using a 
{zfssoip — Yfio5m)AB > 1.0 cut should then efficiently se- 
lect galaxies at z > 6.5. While galaxies at z > 8.0 will 
have colours consistent with (zfssoip — Yfiosw) ab > 1.0 the 
V/ios^-band filter is increasingly contaminated by the break 
reducing the recorded flux and thus the effective selection 
volume. 



3.2 Contamination 

While a single colour criteria is capable of discriminating 
high-redshift objects from the bulk of the galaxy population 
it is still liable to several types of contamination, includ- 
ing both objects with intrinsically red colours and transient 
phenomenon. 

Intrinsically red objects predominantly include galaxies 



1 An approximation to the slope expected for a young star form- 
ing population with no dust 
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Figure 1. The (z/85oz p — 1/io5id)ab colour evolution of a galaxy 
assuming 

fx = \ fi at A rest > 1216A, f x = 0.01 x A" 2 at 912A < 
A rcs t < 1216A (i.e. intervening neutral Hydrogen absorbs ~ 99% 
of the emitted photons) and fx = 0.0 at A res t < 912A. The three 
lines denote /3 = 1.0, —2.0 and —5.0 from top to bottom and 
seek to highlight that the redshift range probed by single colour 
selection is sensitive to the underlying continuum slope. 



Figure 2. Horizontal bars denote the redshift range over which 
each filter probes the rest-frame UV continuum from 1250 — 
3000 A. The vertical bars denote the rough redshift range for each 
dropout sample. 



with prominent Balmer/4000A breaks located at lower red- 
shifts than the target population (for example these galax- 
ies at z ~ 1.5 would contaminate a pure z/gsoip — V/ios™ 
selected sample). Galaxies with prominent Balmer/4000A 
breaks however typically have intrinsically very red spec- 
tra at wavelengths longward of the break location, and can 
be idenitfied using additional colour information. Cool, low- 
mass dwarf stars are an additional contaminant (late M- 
class for the u-drop selected galaxies and LT-class for the i 
and 2-drops) . While in some cases these can also be excluded 
on the basis of additional colour information, for certain fil- 
ter combinations this is not possible and stars can only be 
excluded on morphological grounds. However this is prob- 
lematic especially at low signal-to-noise and may also result 
in the removal of real high-redshift galaxies (see Stanway et 
al. 2008, Wilkins et al. 2011). 

Transient phenomenon capable of contamination in- 
clude stars with high-proper motions or supernovae. Tran- 
sient phenomenon are potential contaminants when there 
is a significant gap between the acquisition of images in 
neighbouring bands. In the context of the observations used 
in this work this occurs between the ACS bviz-band and 
WFC3 YJH-band imaging. An object with an observed 
very-red Zfssoip — Yfio5w colour (i.e. no detection in the 
2/850ip-band) may potentially be a transient phenomenon. 
While the brightest transients can be excluded by examin- 
ing imaging at a previous epoch (assuming it has roughly 
the same spectral coverage) fainter objects are less easy to 
unambiguously identify. 

3.3 Probing the Rest-Frame UV Continuum of 
high-z Galaxies 

To measure the rest-frame UV continuum colours without 
the possibility of contamination from either the Lyman-a 
break or line emission at least two filters above the break 
are necessary - subsequently we refer to these two filters as 
C and D such that \d > Ac- In the context of the avail- 
able WFC3 imaging, the requirement of two filters safely 
above the Lyman-a break limits us to galaxies selected as 



and Zfssoip dropouts (objects with red (A — B) 
colours, where the A filter is the drop-out filter and the 
B filter is the next available filter moving redward); cover- 
ing approximately 4.7 < z < 8.0. While it is also possible 
to identify even higher redshift l/ios™ an d Jfi25w dropout 
galaxies (e.g. Bouwens et al. 2010, Bunker et al. 2010, Mc- 
Clure et al. 2010, Lorenzoni et al. 2011) there is insufficiently 
deep imaging in bands above the rest-frame Lyman-a break 
to measure the UV continuum properties without the possi- 
bility of contamination from Lyman-a emission or the break 
itself. 

In Figure [5] we show the redshift-range over each of the 
three WFC3 filters and the ACS Zfssoip can be used to probe 
the rest-frame UV continuum from 125(0 to 3000 A. Where 
the ranges of two filters overlap they can be used to safely 
probe the UV continuum free from contamination due to 
the Lyman-a break or Lyman-a line emission. At z ~ 7 
(where z-drop selected galaxies will lie) only the J/i25u, and 
Hfi§Q m filters lie safely above the Lyman-a break and it is 
these filters that we adopt to probe the UV continuum slope. 
At z ~ 6 the Yf 105w , Jfvxhw and Hf 160w WFC3 filters all 
lie within the restframe 1250 — 3000A range and can thus 
be used to probe the UV continuum - however to aid in the 
comparison between different redshifts we utilise only the 
Yfio5w and Jfi25 W filters which cover a similar rest-frame 
wavelength range to the Jfi25 W and -H/ieo™ filters at z ~ 7. 
At z > 8 the J/i2s™ filter extends shortward of the rest- 
frame Lyman-a and thus cannot be robustly used to measure 
the continuum colour due to the possibility of contamination 
from the break and/or Lyman-a line emission. 

The redshift range over which individual filters can be 
used to safely probe the rest-frame UV continuum has im- 
portant implications for the specific selection criteria used; 
this is discussed in more detail in the following section. 
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Table 1. Summary of the choice of ABCD niters, where the (A — B) colour probes the Lyman-a break, e.g. for the -z^gsojp-drops 
(A — B) = (z/850Zp — ^/I05u>) an d (C — D) colour is used to probe the UV continuum, redshift range, and the selection parameters 
employed for each drop out sample. "Filter combination used to probe the Lyman-a break. ''Assumes no evolution within the redshift 
range. c The rest-frame central wavelength of the C filter assuming the source is at z = (z). d Number of objects meeting the selection 
criteria in all three fields. 



3.4 Selection Criteria 

While a simple selection criteria based solely on the break 
colour (A — B) can be used to identify high-redshift galaxies 
it is both susceptible to contamination and can introduce a 
number of biases when attempting to study the rest-frame 
UV properties of high-redshift galaxies. Each of these fac- 
tors much be taken into account when designing a selection 
criteria. 

Firstly, contamination from Zouier-redshift objects has 
to be minimised. Figure [3] shows the positions of O — > T 
class dwarf stars (filled points) and the track followed by 
a stellar population formed from through an instantaneous 
burst at z = 10.0 (maximising the strength of the 4000A fea- 
ture at z = 1.5). In order to exclude these objects, while still 
maintaining a reasonable number of objects we introduce a 
cut on the colour used to probe the UV continuum ((C — D), 
see Table [TJ, i.e. (C — D) < pi. The precise value of pi is 
chosen such that for each dropout sample we probe to the 
same maximum value of E(B —V)oaizetu = 0.8. The (C — D) 
colour corresponding to this value of E(B — V) is determined 
by applying the Calzetti et al. (2000) reddening curve to a 
synthetic SED assuming 100 Myr previous constant star for- 
mation, solar metallicity and a Salpeter (1955) IMF. This 
is further elaborated on in Section 4 where we discuss the 
impact of various effects on the UV continuum slope of star 
forming galaxies. The value of pi for each drop-sample (to- 
gether with the other selection criteria parameters) is pre- 
sented in in Table [T] It is important to note that such a 
criterion clearly imposes an explicit bias against very dusty 
galaxies, however as we will see in Section [5] this appears 
to be only problematic if the colour/dust distributions are 
bimodal. To further guard against potential contamination 
we also introduce a criterion that candidate high-z objects 
not have a detection in any filter which probes rest-frame 
wavelengths of < 912A (below the Lyman limit). 

The second consideration is to ensure that the redshift 
range (and mean redshift) probed is not sensitive to the 
absolute luminosity or intrinsic UV continuum colour, and 
further, that for a given absolute luminosity the effective 
volume survey is as uniform as possible. A related consider- 
ation is that the theoretical redshift range should be com- 
patible with the range over which the available filters can 
be used to safely probe the UV continuum. For example, as 
noted in the previous section, it is necessary to restrict the 
z-drop redshift range to z < 8.0 so that the J/i25 W is not 
contaminated by either Lyman-a emission or the break. 



2 1250A is chosen to provide sufficient distance from Lyman-a. 



The first step in this procedure is to determine the re- 
lationship between the break colour (A — B) at a given red- 
shift and the uncontaminated UV continuum slope colour 
(C — D). This is necessary as objects with different rest- 
frame UV colours will, as a function of their absolute mag- 
nitude defined at some arbitrary point on the rest-frame UV 
continuum, enter a single colour selection at a different red- 
shift. To understand this more clearly, it is useful to consider 
the case of two galaxies with the same rest-frame absolute 
magnitude measured at 1500A (M1500) but different spectral 
slopes; one with a red slope (/3 = 0.0, where f\ = A' 3 ) and 
the second with a blue slope (/? = —2.0). If both galaxies 
have a sharp break at A = 1216A, blueward of which the flux 
density is reduced to 1% of that assuming no break, filters 
located to straddle the break would record different colours 
for each galaxy. If candidates were selected solely on the ba- 
sis of the single colour, the redshift range for galaxies with 
one spectral slope would differ from that for another, poten- 
tially resulting in a continuum colour bias. This can be seen 
in Figure Q] where the z/sscup — ^fios™ colour evolution track 
is shown assuming three values of j3 (j3 G {1.0, —2.0, —5.0}). 
Each track crosses a given Zfssoip — Yfiosw value at a differ- 
ent redshift. 

To account for this effect we introduce a criterion that 
relates the break colour (A — B) to the observed colour used 
to probe the continuum (C — D) at a given redshift, i.e. a 
criterion which requires that the (A — B) > f(C — D),. To 
ascertain this relationship we use a range of UV continuum 
slopes and determine (A — B) as a function of the (C — D) 
for a given redshift. This redshift ultimately forms the lower- 
boundary of the selection range and is chosen such that the 
selection window is sufficiently removed from the majority 
of the contaminant population. The relationship between 
(A — B) and (C — D) is found to be approximately linear 
thus we introduce a criteria: (A — B) > P2 x (C — D)+ps. The 
selection window in (A — B) — (C — D) colour space resulting 
from both this criteria and with the (C — D) < pi criteria 
is shown in Figure O The shape of the selection window has 
the additional benefit in that for a given redshift boundary it 
reduces the amount of contamination compared to assuming 
a simple (A — B) based criteria. Indeed, the colour criteria 
often employed to optimise the number of robust candidates 
typically has a similar shape. 

While these two criteria can be used to effectively fix 
the lower limit of the redshift range, the upper limit remains 
sensitive to the UV continuum colour, and for the z/ssoip 
drop galaxies is above the redshift at which the J/1251U filter 
can be used safely. The solution to this problem is to intro- 
duce a further colour criteria, this time using the (B — C) 
and (C-D) colours (i.e. (B - C) < p 4 X (C-D)+p 5 ). The 
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Figure 3. (C — D) — (A — B) colour-colour (where the (A - B) 
colour probes the Lyman-a break, e.g. for the z^g5Qi p -drops 
(A - B) = (zfssoip ~ Y fW5w) and (C - D) colour is used to 
probe the UV continuum, see Table [TJ figures for each of the 
dropout samples (v, i and z from the top respectively) showing 
the selection window (grey), the positions of O — > T class dwarf 
stars, the track taken (over a narrowed redshift range where the 
break colour probes the Balmer/4000A break) by an instanta- 
neous burst of star formation at z = 10, and the location of our 
candidate objects. The size of the symbol is an indication of the 
absolute magnitude M1500 and error bars are omitted for clarity. 



parameters p4 and ps are determined in a similar technique 
to P2 and pz. 

In addition to the colour criterion outlined above we 
also introduce an abolsute magnitude M1500 criteria such 
that M1500 < —18.5. M1500 is determined using the abso- 
lute magnitude inferred from the apparent C magnitude as- 
suming the source lies at the mean redshift and applying a 
correction using the (C — D) colour to scale it to 1500A. 
Given that the pivot wavelength of the C-filter at the mean 
redshift of each drop-out sample is ~ 1500A this correction 
is very small. For this criteria we choose M1500 < —18.5 cor- 
responding to zfssoip = 27.89, Yfwsw = 28.1, J/i25» = 28.5 
at z = 5.0, z = 5.9 and z = 7.2 respectively. For the f/606™ 
and i/775iu selected drops this is significantly brighter than 
the 5<t point source limit, however for the z/ssoip-drop se- 
lection this limit is close to the 5a point source limit sug- 
gesting large uncertainties and incompleteness. Indeed, as 
will be shown in the proceeding section, where we employ a 
series of photometric scatter simulation this appears to be 
the case, and reference is made to this when we analyse our 
results. 

In summary, the selection criteria we employ can be 
defined as: 
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together with a non-detection (< 2a) condition in the rel- 
evant optical filters. These niters, and the selection criteria 
parameters (pi) are given for each dropout sample in Table 

m 

To verify our selection criteria and investigate any fur- 
ther biases we conduct a series of photometric simulations; 
these are described in the following section. 

Finally, it is important to stress that these selection 
criteria are designed to minimise any potential biases in the 
measurement of the UV continuum slope. A result of this is 
that the number of candidates selected using this criteria is 
potentially much smaller than a criteria designed solely to 
maximise the number of high-redshift galaxy candidates. 



3.5 Simulated Redshift Distributions and 
Effective Volumes 

The presence of photometric noise will have the effect of 
scattering objects both in and out of the selection window, 
possibly distorting the redshift distribution and affecting the 
effective volume. Such effects may introduce an additional 
colour bias (other than that explicitly denned in the selec- 
tion criteria, i.e. E(B — V^calzetti < 0.8) or affect the mean 
redshift. 

To investigate this we perform a series of simulations in 
which we populate the original images with synthetic sources 
with various properties (common properties including: IGM 
absorption as given by Madau et al. 1996, and sizes and pro- 
files consistent with the secure candidates from Wilkins et 
al. 2011). Using the same techniques (i.e. use of SExtrac- 
TOR, aperture photometry, and selection criteria) that are 
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used to identify candidate high-z objects we assess the abil- 
ity to recover objects based on their UV colours, redshifts 
and absolute magnitudes. An example of this is shown in 
Figure [4] this figure shows the recoverability, i.e. the frac- 
tion of sources inserted that are recovered using our selec- 
tion criteria, as a function of redshift for various rest-frame 
UV absolute magnitudes (M1500 £ {-21.0, -20.0, -19.0}) 
and UV continuums slopes (/3 € {-2.0, -1.0,0.0, 1.0}). 

For the majority of cases the shape of the recovered 
redshift distribution (and importantly the mean redshift) 
closely matches that which we defined previously, further, 
while the average recovered fraction decreases at fainter lu- 
minosities and redder UV continuum slopes this effect is 
mild with a single exception. In the case of the faintest z- 
drops case we see that both the recovered redshift distri- 
bution is changed (with the mean recovered redshift now 
lower than z — 7.2) and the average recovered fraction is 
decreased by ~ 50 percent. 

The effect of the rest-frame UV absolute magnitude and 
colour in affecting redshift distribution and mean recover- 
ability can be shown more clearly by estimating the effec- 
tive volume for each combination of parameters. To allow 
easier comparison we normalise this to the volume if the en- 
tire theoretical redshift range (as defined by the selection 
criteria) is probed uniformly with a probability of recovery 
of 100 percent. 

The normalised effective volume is shown, as a func- 
tion of intrinsic UV colour, for the three dropout criteria 
and assuming three different absolute magnitudes in Fig- 
ure [5] In the brightest two bins there is only a < 10% 
decrease in the effective volume over the range of colours 
we probe. For the faintest magnitudes considered, and in 
particular the highest-redshift z-drops, the effective volume 
falls significantly (for the z-drops the normalised effective 
volume is: « 0.55 assuming (J/i25 TO — Hfi 60w )AB = 0.0). 
For these faintest magnitudes there is also a significant 
bias against objects with red rest-frame UV colours; for 
example the normalised effective volume of z-drops with 
M 1500 = -19.0 and (J fl2 5m ~ Hf 160w ) A B = 0.6 is w 0.20 
(c.f. (J/12STO - Hj 160w )ab = 0.0: w 0.55). When examining 
the observed distribution of UV colours it is important to 
take these biases into account. 

In Section 15.31 where we consider the observed distri- 
bution of UV continuum colours we also use make use pho- 
tometric scatter simulations in the context of ascertaining 
the scatter in the UV continuum colours we would expect 
purely due to photometric scatter. 



4 THE UV SPECTRAL SLOPE AS A 
DIAGNOSTIC OF STAR FORMING 
GALAXIES 

The UV continuum of star forming galaxies is dominated by 
emission from the most massive, high-temperature (> 3M©, 
OB class) stars. The UV (1216 -> 3000A) spectral flux den- 
sity of the hottest of these stars is well characterised by a 
power law with a blue slope (i.e. if fx oc A' 3 then j3 < 0.0). 
This power-law behaviour principally arises because the UV 
continuum (from 1216 — > 3000A) probes the Rayleigh- Jeans 
tail of the black body spectral energy distribution (though 
it deviates from a pure Rayleigh- Jeans behaviour because 
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Figure 4. The simulated probability of recovering a synthen- 
tic galaxy inserted into the original images (i.e. recoverabil- 
ity) as a function of redshift assuming different properties. 
left panel Assuming M1500 = —21.0 and various values (/3 £ 
{-2.0,-1.0,0.0,1.0}) of the intrinsic UV continuum slope /3. 
right panel Assuming = —2.0 and various values (M1500 S 
{-21, -20, -19}) of the rest-frame UV luminosity. 



of the effects of opacity in the star's atmosphere). At lower 
effective temperatures the range 1216 -+ 3000A probes the 
black body outside the Rayleigh- Jeans tail toward where the 
black body distribution peaks. This effect, combined with 
variation in the opacity (as a function of temperature etc.), 
results in the SED of the cooler, lower-mass stars having a 
redder UV continuum than that of stars with higher tem- 
peratures/masses. 

The sensitivity of the UV continuum of a star to its 
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Figure 5. The relative (to the maximum achievable) effective vol- 
ume probed as a function of the UV continuum slope for different 
absolute magnitudes for each drop sample for the combination of 
HUDF, P12 and P34. 



surface temperature (and opacity) , and thus metallicity and 
mass means that the continuum of a composite stellar popu- 
lation is then sensitive to both the distribution of the masses 
of the stars and the metallicity. The distribution of masses is 
in turn determined by star formation history (SFH) and ini- 
tial mass function (IMF); thus both these factors affect the 
observed UV continuum colours. In addition extrinsic prop- 
erties such as the amount and composition of any gas and 
dust (either in the host galaxy or in intervening systems) can 
produce an effect. In particular the strong wavelength sensi- 
tivity of typical dust reddening curves can cause potentially 
large deviations in the UV continuum colours of galaxies. 

In this section we investigate the magnitude of each 
of these effects. To aid with the comparison of different 
effects we define a default scenario about which we con- 
sider deviations in the amount of dust, the previous SFH, 
IMF and metallicity. This default scenario assumes 100 Myr 
continuous star formation, solar metallicity (Z = 0.02), an 
IMF described by a two-part broken power with a Salpeter 
(1955) high-mass (0.5 < M/M Q < 120) slope (a 2 = -2.35) 
and a shallower low-mass (0.1 < M/Mq < 0.5) slope 
(qi = —1.30), no dust and is constructed using the Pe- 
GASE.2 (Fioc and Rocca-Volmerange 1997, Fioc and Rocca- 
Volmerange 1999) population synthesis model. For a galaxy 
at z = 7.2 (the mean redshift of the z/gsoip-drop sample) 
this scenario suggests a (Jfizsw — Hfieow) colour of —0.037. 



4-0.1 The P parameter 

To compare observations at different redshifts and/or us- 
ing different filter combinations many studies parameterise 
the UV continuum in terms of index of the power law slope 
index /3 (where fx oc A' 3 ), where the default scenario sug- 
gests /Sdefauit — —2.16. However, some studies differ in their 
conversion from the observed colour to /3; in this work we 
convert the (C — D) colour to j3 by determining the value 
of j3 for a given colour assuming the underlying slope is 
perfectly represented by a power law. By this definition 
P = st X (C - D) - 2.0 where si is equal to 6.21, 5.37 
and 4.28 for the v, i and z-drops respectively. 



4.1 Effects Intrinsic to the Stellar Population 

4-1.1 Initial Mass Function 

As noted, the UV continuum of a star is sensitive to its 
initial mass, with more massive (hotter) stars (e.g. O class) 
typically having bluer spectral slopes than the cooler (B and 
A class). Clearly then, any factor affecting the distribution 
of masses potentially can affect the slope of the composite 
stellar population. 

For example, a stellar population containing only stars 
with masses > 20 Mq (and thus containing only O stars) 
will have a bluer UV continuum than a population including 
stars down to 0.1 Mq. This can be seen in the top panel 
of Figure [6] where the J/i25m — Hfieow colour of a stellar 
population at z = 7.2 is shown as a function of the low 
mass cutoff (mi ow ). For a population containing only the 
most massive > 100 Mq stars the predicted Jfi25w — Hf isow 
colour is w -0.15 (/3 = -2.65). 

The relative number of stars can also be affected by as- 
suming an IMF with a different high-mass slope cti (c.f. for 
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a Salpeter IMF o?2 = —2.35). Increasing a?2, i.e. flattening 
the slope, increases the contribution from the most massive 
stars again resulting in a blueward shift in the UV contin- 
uum colour. For example, increasing 02 to — 1 (i.e. a uniform 
stellar mass distribution) results in J/125™ — -H/160™ ~ —0.09 
(/3 — —2.39). The sensitivity of the J/125™ — -ff/160™ colour 
to Q2 can be seen in the second panel of Figure [6] 

4-1-2 Star Formation History 

The strong initial mass dependence of the main-sequence 
lifetimes of stars means that the present day mass distri- 
bution, and thus UV continuum slope, is sensitive to the 
star formation history. After < 10 Myr of continuous star 
formation a population would would still retain most of its 
original high-mass stars. However, after prolonged periods 
of star formation some fraction of the most-massive stars 
will have evolved off the main-sequence. This reduces the 
relative contribution of these stars to the UV continuum 
resulting in a redder slope. In the third panel of Figure [6] 
the (J/125™ — -H/160™) colour (of a population at z = 7.2) 
is shown as a function of the previous duration of con- 
tinuous star formation. As the duration of previous star 
formation increases the UV colour becomes redder. After 
10 Myr (J/125™ — -H/160™) ~ —0.07 while for longer dura- 
tions, > 1 Gyr, (J/125™ - -H/ieo™) > 0.0. 

Star formation histories which are more stochastic, com- 
prised of bursts and troughs, are capable of producing much 
redder UV continuum colours. For example, an instanta- 
neous burst of star formation followed by ~ 50 Myr of qui- 
escence produces (J/125™ — -H/ieo™) ~ 0.07, while still pro- 
ducing a significant UV luminosity. 

4-1.3 Metallicity 

Decreasing the metallicity in stellar cores results in a de- 
crease in the internal opacity allowing energy to be more 
efficiently transported from the core. This allows a star of 
the same mass but lower metallicity to sustain a higher core 
temperature and thus produce more energy. The sensitivity 
of the J/125™ — Hfieow colour to the metallicity is shown in 
the fourth panel of Figure [6]- decreasing the metallicity from 
Z = 0.02 -> 0.004 results in A(J/i 25 ™ - H fl60w ) * -0.08. 
This is an important effect to take into account given obser- 
vations of sub-solar metallicity UV producing stellar popu- 
lations in galaxies at high-redshift (e.g. Pettini et al. 2000). 

4.2 Dust 

By far the largest potential effect on the UV colour of star 
forming galaxies is the presence of large quantities of dust. 
The strong wavelength dependence of typical (e.g. Calzetti 
et al. 2000) reddening curves (where k\ is a monotonically 
decreasing function of A) results in greater extinction in the 
UV (relative to the optical or near-IR) and the reddening 
of UV colours. The reddening of the UV colours then pro- 
vides a potential diagnostic of the magnitude of the dust 
attenuation (Meurer et al. 1999). 

Applying the Calzetti et al. (2000) curve to the de- 
fault scenario inferred SED produces a relationship between 
E(B— V)caizetti an d the UV continuum (C—D) colour which 
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Figure 6. The effect of changes to the lower stellar mass cutoff 
of the IMF (top panel), high-mass slope of the IMF (2nd panel), 
duration of previous continuous star formation (3rd panel) and 
metallicity (bottom panel) to the J/125™ — Hfi60w colour of a 
galaxy at z = 7.2 from the Pegase.2 population synthesis model. 
The vertical and horizonal lines denote the value and colour for 
the default scenario respectively. 
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is approximately linear, i.e. E(B — V) = d\ x (C— D)+d2 (an 
expression similar to that suggested by Meurer et al. 1999). 
Differences between the relationships for different filter com- 
binations arise because in each case the specific filter com- 
bination probes a unique range of the rest-frame UV con- 
tinuum the values of (d 1 ,d 2 ) are (1.22, -0.05), (1.07, -0.04) 
and (0.88, —0.04) for the v, i and z drops respectively. 

Of course, a significant concern, is whether the Calzetti 
et al. (2000) curve is appropriate for high-redshift star form- 
ing galaxies, given that it was empirically derived from local 
starburst galaxies. At high-redshift (esp. z > 6.) there has 
been insufficient time to produce significant amounts of dust 
through the winds of stellar atmospheres suggesting that any 
dust may be dominated by a supernovae component (Dunne 
et al. 2003) and may therefore have a different composition 
and thus reddening curve. However since the precise nature 
and UV-extinction properties of dust in starburst galaxies 
at this epoch is still poorly constrained, we use the Calzetti 
et al. (2000) dust curve as a first approximation. 

Both the effect of dust (though only assuming the 
Calzetti et al. 2000 curve), and the other parameters de- 
scribed in the preceding section are summarised in Figure 

m 

4-2.1 The Interpretation of UV Colours 

ft is important to stress that these relationships (between 
the UV continuum colour and the colour excess) only hold 
for the default scenario defined earlier in this section (and of 
course are only relevant for the Calzetti et al. (2000) redden- 
ing curve). Any deviation in the underlying stellar properties 
(metallicity, SFH, IMF) will have a systematic effect on the 
dust properties derived from the rest-frame UV colour com- 
plicating the interpretation. Thus without additional infor- 
mation (which is scarce at high-redshift) any observed UV 
colour will be unlikely to have a unique interpretation. For 
example, while a J/i25 TO — H/i60ui = —0.037 is consistent 
with the default scenario (i.e. no dust) it also consistent with 
non-zero attenuation together with lower stellar metallicity 
or a shorter duration of previous star formation. 

One exception to this is if extremely blue UV contin- 
uum colours are observed; while the UV continuum slope 
can be made arbitrarily red by adding additional dust, the 
same is not true in the opposite direction (i.e. the UV 
slope can not be made arbitrarily blue). For example, a 
Jfi25w — Hfwow < —0.2 colour (/3 < —2.9) for z-drop se- 
lected galaxy is inconsistent (assuming the accuracy of the 
Pegase.2 model) with Z > 0.0004 metallicity and any po- 
tential star formation history. A secure observation of such 
a colour may then potentially be an indication of an ad- 
ditional effect (see Bouwens et al. 2010a), perhaps due to 
stellar populations with ultra-low metallicity or a different 
IMF. 

While it is not currently possible to independently and 
accurately constrain the metallicity and star formation his- 
tories of high-redshift galaxies this can be investigated in 
a statistical sense (i.e. for large samples of galaxies) using 
galaxy formation simulations. For example, galaxy forma- 
tion simulations can be used can be used to determine the 
star formation and metal enrichment histories of galaxies ex- 
ceeding some UV luminosity threshold. Using these proper- 
ties combined with a population synthesis model can then be 



used to determine the contribution of the SFH and metallic- 
ity to the UV continuum colours. Of course such a technique 
is reliant on the accuracy of galaxy formation models, never- 
theless it provides a first step in constraining the magnitude 
of these effects and thus ascertaining contribution from dust 
to UV continuum colour (Wilkins et al., in prep). 



5 UV CONTINUUM COLOURS OF 
CANDIDATES 

Using data from the deep imaging fields discussed in Section 
2 and the selection criteria laid out in Section 3 we identify v, 
i and z drop candidate high-redshift galaxies with M1500 < 
— 18.5; the location of each galaxy in the (A — B) - (C — D) 
colour plane is shown in Figure together with the selection 
window and potential contaminants. 

Figure [8] displays the (C — D) colour for each candidate 
as a function of the absolute magnitude at 1500A. The ab- 
solute magnitude is derived from the magnitude in the C 
filter corrected using the (C - D) colour to 1500A. The rest 
frame location of the pivot point of the C filter assuming 
the mean redshift are 1416A, 152lA and 1524A for the v, i 
and z-drops respectively; thus the magnitude of this correc- 
tion is typically very small. However, the uncertainty in the 
redshift introduces a systematic uncertainty on the absolute 
magnitude which is dependent on the (C — D) colour (and 
its uncertainties). 

Throughout this section we employ the observed UV 
continuum colour (i.e. (C — D)ab), the f3 parameter and the 
E(B — V) inferred assuming the Calzetti et al. (2000) red- 
dening curve and various assumptions regarding the stellar 
population properties (the default scenario: 100 Myr previ- 
ous continuous star formation, solar Z = 0.02 metallicity 
and a Salpeter IMF combined with the Pegase.2 popula- 
tion synthesis model). 

5.1 Very High Redshift Candidates 

There has recently been significant interest in the literature 
as to the security of many proposed high-redshift candidates. 
Given the limited data and the extremes to which the data 
is pushed this is unsurpring especially when using selection 
criteria aimed at maximising the number of candidates, and 
highlights the need for additional information (such as the 
detection of Lyman-a emission) to securely confirm the iden- 
tity. 

For three reasons the very-high redshift objects consid- 
ered in this work can be considered as more robust than 
Wilkins et al. (2011). Firstly, the criteria employed to un- 
formly select galaxies irrespective of their UVC colour de- 
mands a redder break-colour (A — B) than employed pre- 
viously, reducing the risk of contamination from the main 
galaxy sample due to photometric scatter (however it does 
increase slightly the risk of contamination from intrinsi- 
cally red objects). Secondly we impose a stricter optical 
non-detection criteria compared to Wilkins et al. (2011) 
in the light of other work in this field and concerns over 
the robustness of faint, ambiguous sources, and finally we 
probe a slightly brighter range of galaxies (again reducing 
the contamination due to photometric scatter and spurious 
sources). 
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Figure 7. A summary of factors affecting the UV continuum slope. Line caps denote the maximum deviation that are possible for 
each parameter. For example, Calzetti et al. (2000) dust is unable to produce Jf\25w — Hfi60w < —0.037 without changing any of the 
properties intrinsic to the stellar population from the default values. The horizontal line denotes the J/i25u, — Hf leo-w colour inferred 
from the Pegase.2 SPS model assuming the default scenario: 100 Myr previous continuous star formation, Z = 0.02, a Salpetcr IMF and 
no dust. 



Indeed we see a decline in the number of candidates pre- 
sented in Wilkins et al. (2011) to 9,4 and 5 in the HUDF, 
P12, and P34 fields respectively. This does not necessar- 
ily mean that the remaining candidates from Wilkins et al. 
(2011) are interlopers as many of the candidates are lost due 
to either being outside the new selection criteria (and thus 
redshift range) or fail to meet our brighter magnitude limit. 

Comparing our candidates to other major studies: 
Bouwens et al. (2010c) and McLure et al. (2011), we find 
a high-level of agreement with 17 of our 18 included ob- 
jects being included in one of the two secure candidate lists. 
Specifically in the HUDF field all 10 of our objects are in- 
cluded in both Wilkins et al. (2011) and Bouwens et al. 
(2010c) and 9 are included in McLure et al. (2011). In the 
P34 field 3 of our 5 candidates are included in primary listing 
of Bouwens et al. (2010c), while the additional two are listed 
as potential. 4 of the 5 objects are included in McLure et al. 
(2011), including the 2 potential objects from Bouwens et 
al. (2010c). In the P12 field 2 of our objects are included in 
the primary listing and the third is included as potential in 
Bouwens et al. (2010c) (McLure et al. 2011 does not analyse 
this field). 



5.2 Luminosity Dependence 

The absolute magnitude (or luminosity dependence) of the 
UV colours can be most easily explored by binning the 
galaxies by their absolute magnitude. In Figure [8] the mean 
(C — D) colour (and E(B — V)caizetu value) for 3 absolute 
magnitude bins (-21.5 < M1500 < -20.5, -20.5 < Mi SO o < 
— 19.5 and —19.5 < A/i 50 o < —18.5) for each of the drop-out 
sample is shown. The values of the mean (C — D) colour, 
E(B — V)caizetu and Pc-d together with the standard de- 
viation are also presented in Table [2] 



For both the v and i-drops we find that the mean 
(C — D) colour becomes is redder in the high and intermedi- 
ate luminosity bins compared to the —19.5 < Misoo < —18.5 
bin (while the brightest u-drop bin is bluer than the inter- 
mediate bin it only contains 2 objects and is thus suscepti- 
ble to large uncertainties). This behaviour is consistent with 
the general trend found by other studies (e.g. Bouwens et 
al 2009,2010a) at these and lower redshifts. For the high 
redshift z-drops the trend is less clear, as the intermediate 
magnitude bin has a mean (C — D) colour bluer than both 
the high and low-luminosity bins. Of course, the analysis of 
the 2-drops is affected by the small number of galaxies avail- 
able. If instead 2 bins are used (—21.5 < A/1500 < —20.0 and 
—20.0 < A/1500 < —18.5) there is the indication of a similar 
trend found for the v and i drops. It is important to remem- 
ber that the faintest z-drop bin is potentially biased against 
red (C — D) colours; the implication of this is discussed in 
more detail in the following section where the observed dis- 
tribution of (C — D) colours is analysed. 



5.2.1 Redshift Evolution 

The top panel of Figure [9] shows the mean E(B — V") for 
galaxies with —21.5 < A/1500 < —18.5 in each drop sam- 
ple highlighting the redshift evolution of the rest-frame UV 
colours. A significant redshift evolution is found whereby the 
mean (C — D) colour (and thus the mean value of E(B — V) 
inferred assuming the default scenario) becomes redder at 
lower redshift. However, this is potentially biased as the 
relative effective volumes for the low and high-luminosity 
varies for each dropout sample (this is seen in Figure [5]). 
Specifically the lower-redshift v and i drops have a different 
theoretical luminosity distribution from the z-drops. 

To avoid this bias, and to investigate the luminosity de- 
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Table 2. The number of candidates n, mean, standard error on the mean a/^/n and standard deviation (in brackets) of the (C — D) 
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Figure 8. The distribution of (C — D) colours and the inferred 
E(B — V) values as a function of absolute magnitude (M1500) 
for candidates with M1500 < —18.5. The three panels arc for the 
three dropout samples with z, i and v the top, middle and bottom 
respectively. The large filled points denote the mean (C — D) 
colour for three bins and labels indicate the E(B— V)calzetti and 
(C — D) colour on the left and right respectively. The large error 
bars denote the standard deviation a of the objects around the 
mean while the small black bars denote the standard error on the 
mean (o/y/n). 



pendence of the redshift evolution, we can split each sample 
into the three absolute magnitude bins previously consid- 
ered; while this allows us to tackle the luminosity distribu- 
tion bias, the uncertainties increase due to the small number 
of galaxies in each bin. In the intermediate magnitude bin 
(—20.5 < M1500 < —19.5) there is again evidence for redshift 



evolution such that at lower redshift galaxies typically have 
redder UV continuum slopes. In the faintest bin both the 
i and z-drop selected samples have colours ~ —0.05 while 
the lower redshift u-drops are slightly redder. The brightest 
bin again shows evidence for some variation but is severely 
limited by the small numbers of galaxies. 



5.3 Intrinsic Distribution of UV colours 

The distribution of rest-frame UV colours for each magni- 
tude bin and dropout sample are shown in Figure [lOl With 
the exception of the most luminous bin ( — 21.5 < Misoo < 
—20.5), where there are too few objects (n ~ 1 — 2) to accu- 
rately quantify the scatter, we find that standard deviation 
of the observed colour distributions to be a ~ 0.13 — 0.22. It 
is interesting to consider how much of this scatter is caused 
by an intrinsic distribution in the properties of the galaxies 
(including dust) and not just due to photometric noise. 

To asses this a series of simulations are performed where 
we introduce synthetic sources into the original images, sim- 
ilar to the technique used in Section [33] In this case sources 
are inserted with a single rest-frame UV continuum colour 
corresponding to a synthetic spectrum assuming the default 
parameters denned in Section [4] We then run the same 
source identification/photometry/selection as used to iden- 
tify our candidates. 

These simulated single-model distributions are shown 
in Figure \W\ as solid histograms for each dropout sample 
and each magnitude bin. The standard deviation of the re- 
covered distribution increases at fainter magnitudes (and 
higher redshift for the same absolute magnitude) reflecting 
the expected increase in photometric noise. 

For the faintest magnitude bin ( — 19.5 < M ° < 

v 1500A 

— 18.5) we find that the standard deviation of the observed 
distribution for the v and i-drops is slightly larger than that 
predicted solely from photometric noise, and as noted in Sec- 
tion [5]2] they are slightly redder than the default scenario. 
This suggests that for these galaxies there is some intrin- 
sic distribution of star formation histories, metallicities or 
dust content. For the faintest z-drop bin both the mean and 
standard deviation of the observed distribution are consis- 
tent with the simulated distribution; i.e. the observed dis- 
tribution is consistent with a single intrinsic colour. The 
inferred narrow range of intrinsic properties suggests that 
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Figure 9. The redshift evolution of the mean value of E(B — V) 
and /3 for candidates with —21.5 < M1500 < —18.5 (top panel) 
and for each of the three Absolute magnitude bins considered in 
Scction l5.2l The size of the points in each panel are proportional 
to the log 10 of the number of objects. The open circles in the 
top panel denote the values of E(B — V) inferred assuming the 
mctallicity evolution predicted from galaxy formation simulations 
(see q5.4.1l . 



population. However, changes to the metallicity and star for- 
mation history cannot produce arbitrarily blue intrinsic UV 
continuum slopes, and it is difficult to produce extremely 
red intrinsic slopes while still recording significant UV lu- 
minosity. 

By assuming an alternative scenario with extremely low 
(Z — 0.0004), though non-zero, metallicity, a short duration 
of previous star formation (lOMyr) and a Salpeter IMF, we 
can place upper limits on the amount of dust attenuation. 
For all (i.e. M1500 < —18.5) the v, i and z drops the upper- 
limits, assuming this low metallicity, on the mean value of 
E{B - V) C aizctti are 0.32, 0.20 and 0.06 respectively (c.f. 
E(B - V) C aizetti assuming Z = 0.02: 0.17, 0.06 and -0.06). 



5.4-.1 Metallicity Evolution 

Over the period spanned by z — 7.7 — s- 4.7 the metallicity 
of the UV producing population is expected to increase as 
core collapse supernovae inject enriched material into the 
interstellar medium. As we saw in Section \4. 1 .31 lowering the 
metallicity shifts the intrinsic colour of the stellar population 
blueward. Thus including the effect of an increasing metal- 
licity will then act to slightly flatten the inferred trend of 
E(B — V). To explore this in more detail we take advantage 
of galaxy formation simulations (e.g. Dayal et al. 2009 and 
Dayal, Ferrara & Saro 2010) that predict the evolution of 
the metallicity of the UV producing population relevant for 
our sample. For galaxies selected as v, i and z-drops, with 
M1500 < -18.5 average metallicities of Z/Z G = 0.22, 0.18 
and 0.12 are found respectively (P. Dayal, Priv. Comm.). 
Taking into account these metallicities in the conversion of 
the observed (C — D) colour to the colour excess assuming 
Calzetti et al. (2000) we now find E(B - V) C aizctti =± 0.25, 
0.14 and 0.01 for the v, i and z respectively. These values 
of E(B — V r )caizetti are shown alongside those inferred as- 
suming the default scenario in the top panel of Figure [!J] to 
highlight the increase. 



the effects of the colour bias that was found to affect the 
faintest z-drop bin should not have a significant impact on 
the observed mean colour. 

For the intermediate magnitude bin (—20.5 < M ° < 

v 150IM 

— 19.5) the disparity between the single- model simulation 
and the observed distribution increases; the standard devi- 
ation of the distribution of the observed galaxies is larger 
than that predicted from the simulations (i.e. <r b s > o"sim). 
As noted in Section 15.21 the mean colours of the v and i- 
drops in these bins are also redder than that predicted as- 
suming the default parameters. Again, this suggests that 
these galaxies have a distribution of at least one of the star 
formation history, IMF, metallicity or dust content. 

5.4 The Physical Interpretation of UV 
Continuum Colours 

As illustrated in Section 14.2.11 the interpretation of the UV 
colours of star forming galaxies is ambiguous as the contin- 
uum is affected not only by dust but also by the recent star 
formation history and the metallicity of the UV producing 



5.5 Comparison with Previous Studies 

5.5.1 Bouwens et al. 2010a 

The analysis of Bouwens et al. (2010a), based on HSTWFC3 
observations of the HUDF, finds a significant luminosity de- 
pendence and, of greater interest, extremely blue (/3 ~ —3, 
(Jfi25w —Hfieow) — —0.2) UV colours of the least luminous 
(A/1500 — —19.0) z-drop sources. These are interpreted as 
evidence of an extremely low metallicity population. While 
we can not exclude a luminosity dependence of the UV con- 
tinuum colour for our z-drop sources we do not find signif- 
icant numbers of extremely blue sources at these faint lu- 
minosities. While we observe some galaxies with f3 ~ —3.0, 
the results of our simulations (Section 15. 3[) suggest these 
are consistent with that expected due to photometric scat- 
ter. The principal difference between our analysis and that 
of Bouwens et al. (2010) is the choice of selection criteria; 
while ours is designed to minimise the UV colour bias that 
of Bouwens et al. (2010) is optimised to identify a maximum 
number of robust z ~ 7 star forming galaxy candidates. 
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Figure 10. Observed (lines) and simulated (solid) histograms of galaxies binned by their UV continuum colour ((C — D), sec Table 
[T] for the corresponding C and D filters for each drop sample) for each drop sample (progressing vertically) and absolute magnitude 
bin (horizontally). The simulated histogram is the result of inserting synthetic galaxies into the original science frame before candidate 
selection. These galaxies have synthetic spectra assembled using the Pegase.2 population synthesis model (assuming lOOMyr of previous 
star formation, Salpeter IMF, solar metallicity and no dust - the default scenario outline in Section 4). The deviation from a single observed 
value is then due entirely to photometric scatter and not a distribution of the intrinsic properties. The small vertical lines at the top of 
each panel denote the mean colours for both the simulated (grey) and observations distributions. The standard deviations of both the 
observed (upper label) and simulated distributions (lower label) are also shown in each figure. The shaded region on the right of each 
panel denotes E(B — V)c a lzctti > 0.8 which we do not probe. 



5.5.2 Dunlop et al. 2011 

A similar analysis, based on HST WFC3 observations, of the 
UV colours of galaxies has also recently been undertaken by 
Dunlop et al. (2011) using a photometric redshift based se- 
lection (see McLure et al. 2011) instead of the colour-colour- 
colour selection utilised in this study. Dunlop et al. (2011) 
finds UV colours over the redshift range 2 = 5 — 7 consis- 
tent with /3 ~ -2. (i.e. (C - D) ~ 0.0). While these results 
are consistent with our analysis for the i and z drops we 
find somewhat redder UV continuum slopes for the lower- 
redshift v drops ({/3)\ v ~ — 1.7±0.3), though note our results 
are consistent within the mutual lcr uncertainty intervals. 



in the choice of selection criteria. In addition however there 
are number of other possibilities; firstly the Bouwens et al. 
(2009) study is largely based on different samples, suggest- 
ing the possibility of cosmic variance. Secondly the Bouwens 
et al. (2009) analysis is based on NICMOS observations (as 
opposed to WFC3), raising the possibility of errors due to 
uncertainties on the respective zeropoints. The use of differ- 
ent filters also means we are probing different ranges of the 
UV continuum and finally there remains the possibility that 
one, or both, studies suffer from some degree of contamina- 
tion potentially biasing the mean UV continuum colours in 
one direction. 



5.5.3 Bouwens et al. 2009 

The study of Bouwens et al. (2009) measures the UV contin- 
uum slope over z ~ 2 -> 6 based on HST ACS and NICMOS 
observations in the vicinity of GOODS-South and North. 
While both our measurements, and those of Bouwen et al. 
(2009) are consistent at z ~ 6, we find somewhat redder 
UV continuum slopes at z ~ 5, especially at lower lumi- 
nosities where we find p ~ -1.89 ± 0.24 at Mi 50 o ~ -19.0 
while Bouwens et al. (2009) finds /3 ~ —2.64. The reasons for 
this difference are unclear but again are likely be due to the 
bias against redder UV colours arising from the difference 



6 CONCLUSIONS 

In this work we have used K5TWFC3 imaging of the Hubble 
Ultra Deep Field (HUDF) and two associated parallel deep 
fields to measure the rest-frame UV continuum properties 
of high redshift galaxies (z > 4.5). 

This sample of galaxies was selected using a Lyman- 
break technique and was split into three redshift intervals 
corresponding to objects dropping out in the i>/606™, */775u> 
and 2/g5o; p -filters. This selection criteria employed is tai- 
lored to uniformly select galaxies for a large range of intrin- 
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sic rest-frame UV continuum colours, thereby intending to 
mitigate any non-explicit biases. 

For both the v and z-drop selected samples we find evi- 
dence, as also reported by other studies (e.g. Bouwens et al. 
2009), for a luminosity dependence in the mean rest-frame 
UV continuum colours. Specifically, more luminous objects 
are found to have redder colours. For the z-drop selected 
sample the trend is less clear due to the small numbers of 
galaxies in each luminosity bin. 

In addition we see evidence for evolution in the rest- 
frame UV colours z ~ 4.7 — > 7.7 such that galaxies at lower- 
redshift (over the same luminosity range) are typically red- 
der. While this is in general agreement with other studies 
we observe a stronger degree of evolution. 

By comparing the observed distribution of colours with 
the results of photometric scatter simulations we also in- 
vestigate the intrinsic distribution of UV colours. We find 
the observed distribution of the least-luminous z-drop se- 
lected sample is consistent with a single intrinsic colour 
corresponding to a dust-free star forming population with 
Z ~ 0.02 and a previous constant star formation history 
of 100 Myr. This contradicts the findings of Bouwens et al. 
(2010) who report extremely blue UV [p ~ -3.0) UV con- 
tinuum slopes. A potential reason for this discrepancy is 
use of a selection critieria biased against red UV continuum 
slopes, and assertion supported by Dunlop et al. (2011). At 
lower redshifts, and increasing with luminosity, the observed 
distribution of colours widens such that it is no longer consis- 
tent solely with photometric scatter (<7 b s > o" s i m ) suggest- 
ing that galaxies at these redshifts have a range of intrinsic 
rest-frame UV continuum colours. 

The common interpretation of variation of the rest- 
frame UV continuum colours is that they are due to dif- 
ferences in the amount of dust attenuation. While this is 
likely to be the predominant explanation the rest-frame UV 
continuum colours of stellar populations are also affected by 
the star formation history and metallicity (and initial mass 
function), parameters which are difficult to independently 
constrain. Ascertaining the contribution of these effects is 
challenging from observations, especially at high-redshift. 
Galaxy formation models combined with population syn- 
thesis models may be capable of providing estimates for the 
extent of these effects allowing the determination of the re- 
maining colour deviation due to the effects of dust. 

HST WFC3 NIR programmes in progress (e.g. the 
CANDELS MCT project), together with ground based NIR 
surveys (e.g. Ultra VISTA) will provide further constraints 
on the rest-frame UV properties of the luminous high- 
redshift galaxy population by increasing the sample size and 
covering a wider range in luminosities. However due to the 
limited spectral coverage available at wavelengths greater 
than 1.8fim (i.e. beyond the Hfieow filter) studies of the UV 
properties of z > 8 objects will likely be limited until the 
deployment of the James Webb Space Telescope (JWST). 
JWST will provide not only the increased wavelength cov- 
erage to probe the z > 8 Universe but also the survey effi- 
ciency to discover and precisely characterise large numbers 
of galaxies at lower redshifts {z > 3). 
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